Interplay of ferromagnetism and triplet superconductivity in a Josephson junction 



P. M. R. Brydon^, Boris Kastening^, Dirk K. Morr^ and Dirk Manske^ 

^ Max- Planck- Institut fiir Festkorperforschung, Heisenhergstr. 1, 70569 Stuttgart, Germany 
^ Institut fiir Theoretische Physik, Technische Hochschule Aachen, Physikzentrum, 52056 Aachen, Germany 
^ Department of Physics, University of Illinois at Chicago, Chicago, IL, USA 

(Dated: February 5, 2008) 

In this paper we extend our earlier analysis of the novel Josephson effect in triplet superconductor- 
ferromagnet-triplet superconductor (TFT) junctions [B. Kastening et al, Phys. Rev. Lett. 96, 
047009 (2006)]. In our more general formulation of the TFT junction we allow for potential scattering 
at the barrier and an arbitrary orientation of the ferromagnetic moment. Several new effects are 
found upon the inclusion of these extra terms: for example, we find that a Josephson current can flow 
even when there is vanishing phase difference between the superconducting condensates on either 
side of the barrier. The critical current for a barrier with magnetization parallel to the interface is 
calculated as a function of the junction parameters, and is found to display strong non-analyticities. 
Furthermore, the Josephson current switches first identified in our previous work are found to be 
robust features of the junction, while the unconventional temperature-dependence of the current is 
very sensitive to the extra terms in the barrier Hamiltonian. 

PACS numbers: 74.50.+r, 74.45. +c, 74.78.-w 

I. INTRODUCTION 

A fundamental aspect of the superconducting state is the orbital symmetry of the order parameter, which is often 
regarded as a key indicator of the physical mechanisms underlying the pairing. Although the high-Tc cuprates are 
undoubtedly the best known example of a class of superconductors with unconventional (i.e. not s-wave) order 
parameter symmetry,!^ this is also believed to be a feature of many heavy fermion and organic superconductors. Of 
particular interest is the case of an odd-parity (p-wave, /-wave, etc) orbital pairing state, as this implies that the 
Cooper pair is in a triplet spin state. This of course opens the possibility of exotic magnetic properties of such 
superconductors, p-wave superconductivity was first anticipated shortly after the development of conventional BCS 
theoryjl^l only in the last decade, however, have the first examples of s pin trip let superconductors been discovered. 
The most promising candidates for triplet superconductors are Sr2Ru04,'^'^J^ and UPtsj'^it has also been proposed 
for a number of other compounds, such as (TMTSF)2PF6,1^ UGe2,L5 and URhGe.E^l Unique among these materials is 
Sr2Ru04, as it is well known that its normal state can be described by Fermi liquid theory.SEl 

It has long been known that bound states can form at the surfaces of superconductors or at their interfaces with other 
materials.^^ In a junction where two superconductors are separated by a barrier of sufficiently small width, these surface 
states overlap forming so-called Andreev bound states: this is of particular relevance in the theory of ballistic transport 
through Josephson junctions, as the tunneling through the Andreev bound states dominates the low-te mpera ture 
transport. ^■^ Andreev bound states are also formed in junctions involving unconventional superconductors.l^ESl The 
Andreev states are responsible for many of the unique features of the current through such junctions, due to their 
strong sensitivity to the pairing symmetry of the superconductors on either side of the junction: a well-known 
experimental consequen ce of this sensitivity is the low-temperature anomaly in the Josephson current between two 
d-wave superconductors .I^^EHxhere has recently been much interest in studying Josephson junctions involving p-wave 
superconductors, as the current through the Andreev bound state s is predicted t o have unique characteristics which 
may be considered to be the signature of the p-wave pairing state.^^IISIlSEoElMl xhe use of such superconductors in 
Josephson junctions is therefore expected to produce new phase sensitive devices. 

The construction of novel Josephson junctions also extends to the choice of tunneling barrier between the two 
superconductors. In particular, there has been much interest in the theory of the Josep hson ef fect between two con- 
ventional superconductors separated by complex heterostructures or magnetic materials .'22l2il2ll An excellent example 
of the unusual properties of the latter class of junctions is provided by the predictiorP^l and subsequent experimental 
verification'^ that a sign change in the current as a function of the temperature is possible for a metallic ferromag- 
netic barrier. Reversal of the current across a ferromagnetic bar rier w ith potential scattering is also predicted to 
occur, although the origin of this effect is fundamentally different p£l22l ^jig metallic case, the sign-change is due 
to the temperature-dependence of the decay and oscillation lengths of the superconducting order parameter within 
the barrier; in the case with potential scattering, the effect is produced by the temperature-dependent changes in the 
occupation of the Andreev states. 

Even though they are expected to have an intimate connection, the interplay between magnetism and p-wave 
superconductivity remains poorly understood. A promising route of investigation into this fundamentally interesting 
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problem is the fabrication of devices that combine these two phenomena in a controlled manner I^Sllll in Ref. 
the ballistic tunneling through a Josephson junction constructed by sandwiching a ferromagnet between two Pz- 
wave superconductors was studied, the so-called "triplet superconductor-ferromagnet-triplet superconductor" (TFT) 
junction. The Josephson current (/j) through the TFT junction demonstrated a very rich dependence upon the 
relative orientation of the d-vectors of the two superconductors and the ferromagnetic moment of the tunneling 
barrier, which were all assumed to be parallel to the junction interface. As in the case of s-wave superconductors, 
the ferromagnetic barrier is responsible for a reversal of the current with increasing temperature. The authors also 
pointed out that in certain circumstances the current is very sensitive to the alignment of the junction components, 
and noted that this sensitivity could be exploited to create "current switches" : very small changes in the alignments 
can cause large increases in the magnitude of Ij, effectively tuning the junction between "on" or "off" states, or 
alternatively an abrupt reversal of the direction of current flow can be produced. This has particular relevance to 
the field of quantum information technology, as these "two level"-type effects open the way for the development of 
novel types of quantum bits.f^S Of course, the dependence of the current on the relative orientation of the junction 
components would also act as an important test of the p-wave symmetry of the superconductors. 

In this paper we extend and elaborate upon the work presented in Ref. by the inclusion of additional scattering 
terms in the barrier, in particular a potential scattering term and a coupling to a component of the magnetic moment 
normal to the barrier interface. Our first objective is to assess to what extent the "current switch" effect and the 
sign reversal of the Josephson current with increasing temperature are robust to the more general description of the 
barrier. We show that the "current switch" effect is quite robust, while the temperature-dependent sign reversal of 
the current is considerably more sensitive with respect to the inclusion of other scattering terms. 

Our second aim is to study the emergence of novel phenomena arising from the inclusion of these additional barrier 
terms. Specifically, we predict three new effects at T = 0: (i) if the d-vectors of the left and right superconductors are 
not aligned, it is possible to generate a non-zero Josephson current even for zero phase difference between the two 
superconducting condensates; (ii) for appropriately chosen barrier portentials and angle between the d-vectors, there 
is no net current flowing through the Andreev bound states for a flnite range of phase differences between the two 
condensates; and (iii) the presence of potential scattering terms in a magnetic barrier can substantially enhance the 
Josephson current flowing through it. We furthermore calculate the spin transport properties of the junction in the 
case when spin-flip scattering is absent from the barrier. We flnd that the z-component of the spin current flows even 
for (p — 0, so long as the d-vectors of the two superconductors are not aligned. Finally, we examine the dependence 
of the critical current on the barrier parameters and the alignment of the d-vectors. We present a prediction for the 
"phase diagram" of the TFT junction (as a function of the barrier parameters and the d-vector alignment) in which 
the phases correspond to different locations of the critical current in the current vs. phase relations of the junction, and 
phase boundaries are given by non-analyticities in the critical current. Apart from the inherent theoretical interest, 
this would also make for an excellent experimental test of our knowledge of the TFT junction. 



In the first part of our paper we introduce in detail the Hamiltonian description of the TFT junction [Sec. II A and 
the construction of the associated Bogoliubov-de Gennes (BdG) equations [Sec. II B . After introduction of specific 
interaction terms in Sec. |II C| for the barrier and the bulk superconductors, the BdG equations are solved in Sec. |IID[ 
where we obtain a general expression for the Andreev bound state energies Ea^b f^nd the Josephson current Ij. These 



results form the basis of Sec. Ill where we discuss the dependence of Ea.b and Ij upon the several parameters 



characterizing the system. In Sec. Ill A we focus upon the case of aligned d-vectors, including a discussion of the 
temperature-dependence of the current in Sec. |III A 1[ This is followed in Sec. |IIIB| by the more general case of 
non-aligned d-vectors. The critical current through the Josephson junction, and the dependence of its first-order 
non-analyticities upon the junction parameters, is presented in Sec. IIIC We conclude in Sec. IV with a summary of 
our results and an outlook for further work. 



II. THEORY 



A. General Hamiltonian for a TFT junction 



We derive the BdG equation for the case of a one-dimensional Josephson junction oriented along the z-axis, con- 
structed by sandwiching a ferromagnetic layer of width d between two triplet superconductors [a schematic diagram 
of the TFT junction is provided in Fig. ([iJ]. For the triplet superconductors to the left and right of the junction we 
assume that their respective d vectors d^ and dji lie in the spin x-y plane, and are parameterized by rotation angles 
and with respect to the x-axis. The magnetization of the barrier has a component Mj^ that lies in the x-y 
plane and is parameterized by an angle a with respect to the x-axis, and also a component M|| in the z direction. 
The barrier is also assumed to contain a potential scattering term. 

The Josephson junction is described by the Hamiltonian H = J dz'dzTi.{z' , z), where the Hamiltonian density is 
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FIG. 1: (color online) Schematic diagram of the TFT junction studied in this work. The figure shows the specific choice of 
parameters that we adopt from the end of Sec. |II D] onwards. 



defined by 

n{z\ z) = Ho{z', z) + H±{z', z) + H\\ {z', z) + Ht{z', z). 
The first term in Eq. ([T]) describes the kinetic energy and the potential scattering by Lip{z), 
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^l + Up{z) tpcriz), 



(1) 



(2) 



where is the annihilation field operator for an electron of spin a. In general, we can have different effective 

masses m on either side of the junction and within the barrier itself, respectively m^, mp and ms- This situation 
may be included in Eq. Q by assuming a z-dependent m. 

The interaction between the magnetic moment of the barrier and the spin of the conduction electrons is given in 
the standard way via 



Wi(z',z) + 7^,[(z',z) = -M(z',z) •^^/.i(z')o•a;3^/'/3W• 
Q,/3 



(3) 



The term Ti.± describes the scattering of the quasiparticles by the magnetic moment in the x-y plane. The orientation 
of the moment with respect to the a;-axis is given by a, and so we have the vectorial representation for the scattering 
potential M±{z',z) = S{z — z')A^j^(z)(cos a, sin a, 0). We thus obtain 



n^{z\ z) = ~5{z - z')M^{z) e-^>|(z')^i(^) - e^"VT(^')V'I(^) 



(4) 



The scattering of the quasiparticle by the magnetization in the z-direction is described by the potential Afj|(z',z) = 
S{z — z')A4\\ (^)(0i 0, 1), and so we have 



H||(z',z) = -^(z-z')-M,i(z) [v^|(z')^TW-V'I(^')V'iW 
Finally, we have the triplet pairing term 

Hriz^z) = A{z',z) ( ^l(z') ) (d - ^)Ha2) ( ^jj^j j + h.c. 



(5) 



(6) 



with complex d- vector d — {d.j;,dy,dz). We choose real d- vectors (so that they are "unitary", i.e. d x d* = 0) that 
lie in the x-y plane, d = (cos 0, sin 0, 0). For this choice of d- vector, the pairing term may be re-written as 



nT{z',z) = {A(z',z)[e-'V|(z')^|W-e^VI(^')V'IW 

-A* (z', z) [e^ Vt (^') V^T (^) - Vi i^'Hl i^)] } 



(7) 
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where we have used the relationship A(z', z) = —A{z, z') to obtain the last line. 

For performing a Bogoliubov transformation it is convenient to introduce a matrix notation. We define the spinor 



vl/(z') = ( ^|(^') ^t(^') M^')} 
which obeys the matrix anti-commutation relation 

{^{z),¥{z')} ^ S{z- z')i. 
We may then write the Hamiltonian in the form 



H ^ 



1 



dz'dz^^z')n{z',z)^{z). 



The matrix Hamiltonian density Ti.{z' , z) is defined by 

7i:(z',z) 



f 6,,, {T+Up-M\\-fi) e-'^A{z\z) -5,,.M±e-"' 

-e''>A*{z',z) -5^^. {T+Up-M\\-n) 

-6zz'MA_e'°' 5^^, {T+Up+M\\-n) 

^33'Al±e-'" e-^'>A*{z\z) 



(8) 



(9) 



(10) 



\ 

-e*«A(z',z) 

-5,,. (r+z^p+M||-A*) y 



(11) 



where we have adopted the abbreviations Szz' = S{z — z'), T = — Up =Up{z), A4± = Ai±{z) and = Ai\\{z). 
We can easily see from Eq. (11) the justification for the separate parameterizations for the longitudinal and transverse 
components of the barrier magnetization: the former do not mix the spin states of the quasiparticles, whereas for 
the latter the z-component of spin is no longer a good quantum number. This has important consequences for the 
Andreev states and the Josephson current. 



B. Bogoliubov transformation and Bogoliubov-de Gennes equations 



To diagonalize the matrix Eq. (Ill we perform the unitary transformation 

n 

= J dzAtiz)^{z), 

where 

{all O-n Pi, Pn )^ , 



and 



I <,„(2) M/3,n(^) V*p^^{z) \ 

Va,n{z) KA^) Vfj^n{z) u'p^^{z) 

Wc,n{z) w/3,n(^) a;* „(z) 

{z) <,„(2) xp^^{z) w;^^„(z) / 



(12) 
(13) 

(14) 
(15a) 



and A+(z) is the Hermitian conjugate of An{z). The diagonalizing spinors ri„ obey the matrix anti-commutation 
relations 

^ 5„irii. (16) 

From the anti-commutation relations Eqs. ([9]) and we obtain the orthogonality relations for the matrix Am{z) 

dzi+(z)i„(z) = j dzA.m{z)A+{z) = SrarX, (17) 

^i„(z)i+(z') = Y.^n{z)An{z') = 5{z-z')i. (18) 
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Using Eq. (12 1 we re- write the Hamiltonian in terms of the ri„, 

H=^j dz'dz¥{z')n{z',z)^{z) ^^J2^L J dz'dzA+{z')n{z',z)An{ 

m.n '- 

Since the Am{z) are assumed to diagonahze 'H{z' , z), we require that 

dz'dzA+{z')H{z',z)A^{z) = EJnm, 



nr. 



(19) 



(20) 



where En is a diagonal real matrix. Using the orthogonality condition Eq. (18 1, after some algebra we obtain the 
most general form of the BdG equation^^U 



dz'n{z,z')An{z')^ An{z)En. 



(21) 



From the definition of A„(z), we may recast this expression as the more familiar eigenvalue problem. Defining the 
vector of c-numbers 



*j\n(2:) = ( Ui.n(z) V3,n{^) Wj^n{z) Xj^n{z))'^, 



where j — a, P, we may write the first (j — a) and third (j = (3) column of Eq. (21 ) as 

rfz'7i(z',z)$,-„(z') = £;,■„$,■„(/). 



(22) 



(23) 



The second and fourth column of Eq. (21 1 are obtained from these by complex conjugation. We suppress the a and (3 



subscripts in Eq. ( 23 ) , as the equations for j — a and j — P are identical. Since we are interested only in the energies 



of the Andreev bound states (i.e n = 0), we also drop the subscript n, assuming henceforth that we are referring only 
to the n = states. That is, we make the replacements Ejn E, <I>j„(z) ^{z) in Eq. (23). 



C. Specific interaction terms 



Up to now our discussion has been very general: in Eq. (11) the forms of the scattering potentials and the barrier 



width are left indeterminate. For the remainder of this paper, however, we consider only the case where the barrier be- 
tween the left- and right-superconductors is a point contact, and so we adopt the physically reasonable approximation 
that it is of infinitesimal width.ESl We therefore replace the scattering potentials by (5-functions: 



Up{z) = UpS{z), M±{z) = Mj_S{z), 7W|i(z) = A/,i(5(z) 



(24) 



Without loss of generality we choose M± > 0. 

For concreteness, we assume a gap with p^-symmetry. In momentum space the superconducting gap therefore has 
the form 



A(fc) = At sin ka, 



(25) 



where a is the lattice constant of the system. It follows from Eq. ( 25 1 that in real space the gap takes the form 

At 



A{z',z) 



2ia 



[S(z' — z + a) — S(z' — z — a)]. 



(26) 



Consider now a function f{z) that is slowly varying on the length scale of the lattice. We may then approximate the 
integral as 



dz'A{z', z)f{z') = -^[f{z -t- a) - f{z - a)] « -zAt9./(z). 



(27) 



This is relevant for our situation, as the components of <i>(z) are assumed to vary only on the scale of the supercon- 
ducting coherence length. A^ is different on the left (z < 0) and right (z > 0) sides of the junction and given by 
A^e**^ and Afle*"^^ respectively, with Al._r real and positive. 
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D. Solution of the Bogoliubov-de Gennes equations 



Under our assumption of (5-function potentials, and using the approximation Eq. (27 1 for the superconducting gap 
in real space, we may simplify the BdG equation Eq. (23 1 to 



where 



H 



( T+{Up-M\\)5{z)-^ 

-e'"Afj_(5(z) 
\ 



-ie-^'^ATd,. -e~'°'Mj_S{z) 
-T-{Up-M\\)5{z)+n 

T+{Up+M{)5{z)-^x 





e'"Mj.(5(z) 
-T-{Up+M{)5{z)+^i ) 



(28) 



(29) 



and T is defined as for Eq. (111. The BdG wavefunction <i>(z) obeys two boundary conditions at the interface. If 
<i>^(z) is a solution of Eq. ( 28 ) in the v — L^K superconductor, the first boundary condition is given by the continuity 
of the wave-function at the junction. 



*l(0) - $fl(0). 



(30) 



In order to obtain the second boundary condition, we integrate the BdG equations across the junction barrier at z = 
by applying the operator J^^* dz and subsequently letting ^ — > 0. The only contributions come from singular terms, 
which are the second derivatives of the components of and the (5 functions modeling the barrier itself. We hence 
obtain the condition 



9.*i?(o)-a.$L(o) 



{Z — g')ao — g{ao cos a — ia^ sin a) 

-(/((jQ cos a + i(T3 sin a) (Z + g')aQ 



1>i?(0), 



(31) 



where is the Fermi momentum on the z^-side of the junction, defined in terms of the chemical potential /x^ 
and the dimensionless couplings Z, g and g' are defined by 



Z 

9 
I 

9 



rriBUf 



m,BM\\ 



(32a) 
(32b) 
(32c) 



Note the appearance of the effective mass ms within the barrier region in Eq. ( 32 1 . Although the barrier region is 
approximated to be infinitesimally small, in any realistic situation it will be sufficiently thick to define an effective 
mass. 

Since we are interested in the Andreev bound states, our solution for <&;y(2:) must vanish in the limit \z\ — 00. We 
hence adopt the ansatz 



where 



( Uu,± wi/,± Wu,± Xv^± y 



(33) 



(34) 



are real vectors describing right-moving ^1,.+ and left-moving ^i,.- solutions, Ki, and are real and positive, and we 
have ai, = -|-1(— 1) for v = L{R). 

In the bulk superconductors, we find after inserting Eq. (33) into Eq. (28 1 that the BdG equations simplify to a 
pair of 2 X 2 eigensystems: 



^ 

2m „ 



{ai,K^±ipj,)'^ - ^i^ -ie '^"'> Ai,{cr„Ki,±ipi,) 



2m„ 



{ai,K^±ip„)^ + 



Xu,± 



E 



= E 



(35) 
(36) 
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Both Eqs. (35 1 and (36) yield the same eigenvalue equations 



2mi, 



(37) 



The right hand side of Eq. (37 1 must be real and positive. As the Andreev bound state wavefunctions are exponentially 
decaying in the bulk superconductors, we require that they be subgap solutions, i.e. have energy \E\ < k^A^. This 
is only possible if and are non-zero; by setting the imaginary part of Eq. ( 37 1 to zero we obtain the relation 



A 2 



and therefore 



nil, 



■my ^2 _ ^^^^ 



(38) 



(39) 



This expression may be considerably simplified by observing that in physically realistic situations the maximum gap 
magnitude is much smaller than the Fermi energy, i.e. fc^A^ <C /i,/ = h?'kH2my. Fur the rmore, from the requirement 
that E^ be positive we deduce that <^ki,, hence allowing us to approximate Eq. (39 1 by 



E' 



(40) 



Since we have E"^ < k^Af,, this is a subgap solution, as we require. Furthermore, using the inequalities 
Ky,2myAy/h'^ ^ ky, we see from Eq. (38 1 that « k^; in this limit the solutions of Eqs. (351 and (36 1 become 



ETcTy 



±e-*(«--'^-)fc„A, 



E^a^ 



Note that k^, is explicitly dependent upon E, and we obtain by solving Eq. (40 1 

where vp^ — hk^/m^ is the Fermi velocity on the p side of the barrier. 
From Eq. (40 1 we may define the parameter 7,^ S [0,7r] by 



(41) 
(42) 

(43) 



E 



kiiAu 



sm7y 



This allows us re- write Eqs. (411 and (42 1 in the simple form 



(44) 



(45) 
(46) 



Since only the phase difference between the left and right superconductors is of physical significance, we set 0l = 
and (pR — (j) without loss of generality. 

Substituting Eqs. (451 and (46 1 into Eq. (33 1, and then substituting this into Eqs. (30 1 and (311 we obtain explicit 
conditions on the u^,± and w^ ± imposed by the boundary conditions. These can be summarized in matrix notation: 
defining the vector ^l^uw by 



( UL+ UL~ WL+ WL- UR+ Ur- Wr+ Wr^ )^ 



we then have 



0, 



(47) 
(48) 
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where 



AfxFT 



Mil Mi2 
M21 M22 



contains all the information from the boundary conditions: the top two matrix components refer to the continuity 
condition Eq. ( 30 1 , while the bottom two components fix the condition on the derivative Eq. ( |31| . Explicit expressions 
for the 4x4 matrices Mn, M12, M21 and M22 can be found in App. [A] It is clear from Eq. (48 1 that non-trivial 
solutions for the components of ^{z) can only be found when ilixFT is not invertible, i.e. detMrPT = 0. The 
determinant may be explicitly evaluated as 



(e^"^/64) detMxFT 



- [3 + 3 cos cos 2^b. — cos 27^ — cos 27i^ + 4 cos 20 
8 

+4cos2(^i — Or) + 16cos(^L — da) sin7iSin7ii;C0S(^] 
—%Zg cos 7l cos 7^^ sin(^i — sin 
+ (r + r-i) {2(^2 + c,'2) cos^ 7l cos^ 7^ 

- cos 7l cos 7fi[cos(^L - ^fl) cos (/) + sin 7^ sin 7^] } 

-|-^(r^ + r"^) cos^ 7l cos^ 77? + 4(Z^ - g'^)^ cos^ 7l cos^ 7_r 

— 4(Z^ + 5'^) COS7L cos7fl;[cos(^i, — ^fl) cos (/) + sin7i sin7/j] 
-252(4(^2 - cos^ 7i cos^ 7« 

—2 cos 71, cos7fl[cos(^L + ^fl) COS0 — sin7i sin7flj 

-l-rcos^ 7fl(sin2 7l + cos2^i) + r^^ cos^ 7L(sin2 _j_ cos20fl)} 

+45^ cos^ 7l cos^ 77?, = 0, 



(49) 



where we have adopted the notations r = ki^/kji, 9]^ = 9]^ — a and 9fj — 9ii — a. The energies E' of the Andreev 
states follow from the solution of Eq. (49l together with the definitions of [Eq. (44l] and k^, [Eq. (43l]. As such, 
it must contain the underlying symmetries of the Hamiltonian. In particular, we note that the E depend on 9^, Or, 
and a only through the differences 9l — 9^ — a and 9f( = 9ii — a, as dictated by rotation invariance around the 
z-axis. Without loss of generality, we therefore set 6*^ = and 9r = 9, i.e. the direction of d^^ defines the x-axis. 
Furthermore, we see that the solutions of Eq. (49) are invariant under (5,0;) ± tt), as expected from the 

definition of 7i_L in Eq 
E 



Eq.(49 1 also provides a general condition for the existence of zero-energy solutions. When 
we have from Eq. (4M the identity cos7i^ = COS7L = 0; the determinant then reduces to 



detMxFT = 2cosi(6'- 



l cos : 



= 0. 



(50) 



Thus there is a zero-energy solution whenever cos ^(0 ± 0) = 0. This condition is independent of the value of r and 
the details of the scattering potential at the barrier. 

Although Eq. (49 1 is a general result, the solution for the Andreev states can be considerably simplified by adopting 
the assumptions fc^ = kji, mj^ — m^j and kj^A]^ — k^Aji = kpA^. This is the natural situation when the supercon- 
ductors on either side of the gap are made from the same material. The determinant then reduces to a quadratic 
equation in E^: 



D^k%A^^ 



4A 



E^ 



45^ = 



(51) 



where 



i {(1 + 2.g2 + c/'2 + Z^) + (1 + + Z^) cos ^!)COS 9 
cos{9 ~ 2q!)[cos0 - cos 0] + 2Zg' sin sin 0}, 



B 
D 



icosi( 



) COS i ( 



v/(l + .g2 + .g'2 -Z2)2 



4Z2 



The positive solutions of Eq. (51 1 define the Andreev bound states: 



E, 



kpAoVD 



Vda + b + {-)Vda - B 



(52) 
(53) 

(54) 



(55) 
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We consider Eq. (55 1 to be the central result of our analytic consideration. As shown in App.jB] this expression can 
be further simplified for the case treated in Ref. \22\ in this situation it is also demonstrated that real solutions exist 
for all parameter values. 

The last remaining step is the calculation of the Joscphson current //. This is given bjff2] 

^J^-nT.^ t^MEi/2kBT), (56) 

l=a,b ^ 

where is Boltzmann's constant. Since we are most interested in the low-temperature current, where the transport 
through the Andreev states dominates, using Eq. (55 1 we may obtain an explicit expression for // at T = 0: 



efdEg dEk\^ ekpAp / D djPA+m 

h\d^ d^J h y DA+\B\ del) ' ^ ' 



where, for later reference, we have explicitly 



'^'^ ^l^^ "^^'^^ - 2a) - (1 + + Z2) cos e] ~ sgnS} sin0 

+ ^DZ^'sin^cos^. (58) 



Note that we henceforth work in units where h~\. 



III. RESULTS 



In this section we study the TFT junction in the special case fc^ — kfj, rriL — mn and k^A^ — knAn = kpAo. We 
investigate the dependence of the Andreev state energies [Eq. (55l] and the current through them [Eq. (56 1] on the 
different junction parameters. In Sec. |III A] we discuss the Andreev state energies and current for the case of aligned 
d-vectors, with a discussion of the temperature dependence of Ij given in Sec. Ill A 1 The more general scenario with 
non-aligned d-vectors is presented in Sec. IIIB Lastly, in Sec. Ill C we consider the critical current as a function of 
g, Z, 9 and a for a junction with g' — 0. Unless otherwise stated, all results are for T = 0. 



A. Aligned d-vectors 

We first consider the case when the d-vectors of the left and right superconductors are aligned, i.e. 9 = 0. As was 
shown in Ref. .22. for a barrier with only .g ^ 0, the (/)-dependence of the current depends crucially upon the alignment 
of the transverse component of the barrier magnetization with the d-vectors. For M _L d^^ji, i.e. 

a = Qfc = (2n -|- l)7r/2, n G Z, 

it was found that the Andreev states are degenerate and Eq. ( [55| simplifies to 

Ea,b = fcFAo\/^cos(0/2), (59) 

which is plotted as the black solid line in Fig. (pi (a). This result is identical to the Andreev state energies obtained 
at a potential scattering barrier (i.e. Z ^ only),^° as well as for a barrier with only a longitudinal moment {g' ^ 
only), see Fig. ([4|. Because the zero-crossings at 

(j)= (j}zc = C^n+l)?:, neZ, 

occur with dEa.b/dcj) ^ 0, we hence find discontinuous jumps in /,/ [shown in Fig. ([2])(b)]. The degeneracy of the 
Andreev states is lifted when there is a component of the magnetic moment parallel to the d-vectors (i.e. a ^ ckc)- 
This holds for all (f) except for the level crossings located at 

4> — 4>LC = 2ri7r, n G Z, 



as shown in Fig. (|2|(c). At (/) = (j)zc both the a and b states have a stationary point with respect to (f>, and the 
discontinuity in the current is therefore removed [Fig. (p])(d)]. 
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(c) (d) 

FIG. 2: (color online) ^-dependence of (a) the Andreev states and (b) the corresponding Josephson current for g = 1.0, 6 = 
and a = 7r/2; t/i-dependence of (c) the Andreev states and (d) the corresponding Josephson current for g = 1.0, 6 = and 
a = tt/4. 



Including a finite potential barrier Z when g ^ has a very significant effect upon the (/)-dependence of the Andreev 
states, in particular it lifts the degeneracy of the states ai a — ac and removes the level crossings at = (jj^c for 
a ^ ac [Fig J2])(a) and Fig. ^{c) respectively]. In the latter case, the ^-dependence of Ea is strongly modified by 
Z: in Fig. (|2|(c) we see that not only does Ea exhibit almost sinusoidal variation with ^ at Z = 1, but also the 
sign of dEa/d(j) is reversed relative to the Z = case in the vicinity oi 4> = 4>lc- This has interesting implications 
for the Josephson current, which displays a moderate enhancement above the Z = Q values for all 0, as can be seen 
in Fig. ([2|(d). Because of the modification of Ea near 0lci the contribution to Ij from the a and b states is of the 
same sign for all values of 0. Furthermore, although \Ei,\ and hence also the current contributed by the 6 state is 
reduced by a finite Z [Fig. ([3])(a)], the increase in \Ea \ for Z < g can compensate, leading to the over-all enhancement 
of I J at all (p values. At higher values of Z, the monotonic decrease in |£'q| causes a reduction in Ij and we therefore 
obtain the current maximum shown in Fig. (|3])(b). In contrast, for a — ac the maximum in Ij is very much reduced, 
as there is no change in the curvature of the Andreev states and so the maximum is due solely to the maximum in 
\Ea\ (not shown). 

The change in the 0-dependence of the Andreev states with the inclusion of a finite g' is much less dramatic than 
for a finite Z. The degeneracy of the Andreev states at a = ac and the level crossings at a ^ ac both remain intact 
when g' ^ 0, as seen in Fig. (|2])(a) and Fig. (|2])(c) respectively. At a = ac the only effect of g' is to narrow the allowed 
range of Andreev state energies, thus leading to a suppression of Ij as seen in Fig. (2])(b). For a ^ Uc the reduction 
in the |-Ea,&| is not uniform, as evidenced in Fig. ([2])(c) by the greater decrease in the magnitude of Ea &t (j) = 4>zc 
compared to that at = 0lc- Although the modification of Ea can lead to very weak increases in // for close to 
4>zci only a monotonic decrease in the current with increasing g' is observed for other values of 4> [e.g. see Fig. (|3|(b)]. 



We now consider the case g = Q. As can be seen from the BdG equation Eq. (28 1, the absence of a transverse 
component of the barrier magnetization means that spin is a good quantum number. This has the interesting 
consequence that the effective barrier potentials can take different values in the two spin channels: for the spin-cr 
electrons there is a combined charge and magnetic barrier of value Z — ag' where a = +(— ) for spin-| (|). We may 
therefore label the Andreev state energies by cr, obtaining from Eq. (Issl) the expressions 



E^ = fcj.AoVA^cos(0/2) (60) 

where = [1 + {Z — ag')^]^^. We thus require that both potential and magnetic terms be present in the barrier 
to achieve a spin-splitting of the Andreev states. This is shown in Fig. (|4])(a), where the Z — states (black lines) 
are split into the states E^ (inner red broken lines) and E^ (outer red broken lines) . Due to the zero crossings of the 
Andreev states, we find discontinuous jumps in the current [Fig. (|4|(b)]; as can also be seen in Fig. (|4])(b), a finite Z 
can give a slight enhancement of Ij when g' ^ 0. This is due to the maximum in | which occurs when the potential 
and magnetic terms in the spin-l sector cancel each other at Z = g' . The current flowing through the junction is 
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FIG. 3: (color online) (a) Andreev state energies and (b) Josephson current as a function of i^, where p is either Z or g' . Note 
that Z — if = g' and vice versa. We assume — 0, (j) = 2tt/5 and a — 7r/4. 




FIG. 4: (color online) (/!>-dependence of (a) the Andreev states and (b) the corresponding Josephson current for g' = 1.0 and 
9 — 0. a is undefined when g — 0; for g ^ 0, we take a = n/4. 



spin-polarized, since the spin-t states contribute more to the current than the spin-| states. As the current flowing 



through each spin sector is a constant multiple 
current polarization P, defined 



/_D|) of the total current, it is possible to speak of a 



JT 



^l + (Z + g'r~^l + iZ-g')^ 



9') 



(61) 



Z^, shown in Fig. (gKa). This maximum 
oo and hence also 



As a function of g' , the polarization takes a maximum at g'^ax 

value of the polarization asymptotically approaches 1 (i.e. fully spin-polarized current) as Z 
5™a. ^oo [Fig. (HKb)]. . . . _ 

When all barrier parameters are non-zero [the blue dot-dashed line in Fig. (E|(a)], the ^-dependence of the Andreev 
states most closely resembles that of the Z ^ 0, = case [Fig. (|2|(c)" Referring to Fig. (|4|(a), we see that 
for a ^ ac not only are the a and b states non-degenerate at ipzci but they also have stationary points there, 
hence removing the discontinuities in /j. At a = a^j the states are degenerate at </) = (j)zc and show simple cosine 
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FIG. 5: (color online) (a) Spin-polarization P of the Josephson current as a function of g for various values of Z. We have 
g = a = = in all lines, (b) Maximum spin-polarization of the current with respect to g as a function of Z. For given Z, 
the maximum in P occurs at g' = g'^ax = Vl + Z^. 
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FIG. 6: Comparison of the dependence of the current on the barrier potentials for 6 = Q, a = n / A and <j) = 2n/5. (a) Current 
as a function of g for different values of g' = Z; (b) Current as a function of Z for different values of g' and fixed g — 1; (c) 
Current as a function of g' for different values of Z and fixed <? = 1. 



dependence on (j) as in Eq. (59 1, but with distinct amplitudes -D6(a) = [I + {Z + {—)\/g'^ + g''^)^]~^ (not shown). 

As displayed in Fig. Q (a), the presence of the other barrier potentials substantially enhances the current for g > 1. 
At smaller values of g, there is only a monotonic depression of the current with increasing g' and Z. We do not 
observe a maximum in Ij as a function of g for any choice of the other barrier potentials. In contrast, increasing Z 
at fixed g and g' [Fig. ([6])(b)] leads to a clear maximum at Z ^ g + g\ with a greater peak current than at g' — 0. 
Defining the maximum current enhancement by the potential scattering Z as 

Szlj = niax{//} - Ij{Z = 0) 

where max^j/j} is the maximum value of Ij with respect to Z with (j) and the other barrier potentials fixed, we 
therefore see in Fig. ^{h) that 5z Ij increases with increasing g'. When Z = 0, the current as a function of g' 
shows only a weak maximum at 5 > 1 which is most pronounced at fixed ~ (j)zc (not shown) ; in the presence of a 
potential scattering term Z > 1, however, \Ij\ develops a clear maximum at g' ^ Z for all (j) as shown in Fig. (|6|(c). 
The current enhancement when all three barrier parameters are non-zero is most pronounced when g' ^ Z ^ g,l; 
in this limit we approach the case g = 0, when the maxima is due to the cancellation of the magnetic and potential 
scattering terms in the spin-f channel, and there is consequently very little g-dependence of the current [Fig. ([6]) (a)]. 
Thus, the maxima in Fig. ([6|(b) and (c) occur in spite of and not because oi g ^0. 

It is interesting to compare the dependence of Ij on g' and Z for g = to the results for a singlet superconductor- 
ferromagnet-singlet superconductor (SFS) junctionP^ although a maximum in the current through the SFS junction 
is found as a function of Z at constant g' , the value of /,/ at this maximum (Z ^ g') vanishes as Z,g' 00. In 
the case here, in contrast, the current at the maximum approaches a finite value in this limit. The difference can be 
understood as arising from the fact that in the singlet case, the superconducting correlations are perturbed no matter 
in which spin channel the scattering occurs. Even if scattering occurs only in one of the spin-channels (as is the 
case for Z — g'), an effective scattering is induced in the other spin channel through the anomalous superconducting 
correlations. As a result, for Z = g' and Z^g' — )■ 00, the tunneling of electrons and hence the Josephson current 
are completely suppressed. In contrast, in the triplet case for d^jj _L M and Z = 17', only the superconducting 
correlations in the spin-| channel (where the effective barrier potential is Z -f 5') are suppressed, while those in the 
spin-l channel are not (here the effective barrier potential vanishes). As a result, the Josephson current remains finite 
in the limit Z, 5' — > 00 and is solely carried by the spin-l Andreev state. 

The current displays a periodic modulation as the magnetic moment is rotated about the z-axis as shown in Fig. ([7|. 
The amplitude of the modulation increases as </) approaches 0zCi with Ij becoming sharply peaked at a — ac for 
(j) (pzc- In the limit — = (f>zc i 0+, the current vanishes unless M _L d^jj. This behaviour can be 
understood as follows: for all a ^ ac, Ea^b has a stationary point at (p — (pzc [see Fig. (|2])(c)], and so Ij vanishes as 
4> — > (f>zc': for ct = ac, however, there is no stationary point at 4> = (pzc and so hm^^^± dEa^b/94' 0; which thus 
gives a finite Josephson current. Close to (p = (pzc, therefore, the current shows a "switch"-like dependence upon 
the orientation of M in the x-y plane. That is, small variations in a can lead to large changes in the magnitude of 
Ij, "switching" the junction from the "on" -state (/,/ 7^ 0) to an "off" -state (/,/ ~ 0). This Josephson current switch 
survives in the presence of a finite Z or g' , although there are significant differences between the two cases. For 
g' 7^ 0, we observe at a = ckc a reduction in the magnitude of both Ij and the amplitude of the current oscillations. 
For Z 7^ 0, in contrast, there is an increase of the Josephson current at all values of a. The enhancement is almost 
a-independent for (p ~ (pLC, whereas for (p ~ (pzc the enhancement is concentrated at the a = ac peak, strengthening 
the switch effect. 
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FIG. 7: (color online) Comparison of the Z- and jr'-dependence of the current as a function of a for = 0. (a) Z- and 
^-dependence of the current for g — 0.5 and g' — 0; (b) Z- and 0-dependence of the current for g — 1.5 and g' — 0; (c) g'- and 
(/>-dependence of the current for g = 0.5 and Z — 0; (b) g'- and (/)-dependence of the current for g — 1.5 and Z = 0. 




(C) T/T^ (d) fTi 

FIG. 8: (color online) (a) Temperature-dependence of the current for various values of g, with a — 0, 9 — 0, (j> = 7i"/2 and 
Z = g' = 0. (b) Temperature-dependence of the current for various values of a. g = 2.0 and other parameters are as in (a), (c) 
Effect of non-zero g' or Z on temperature-dependence of current at a = 0.157r and other parameters as for (b). (d) Andreev 
states as a function of for the parameters in (c). The dotted vertical line corresponds to the situation in (c). 

1. Temperature Dependence 

In Ref. I22l it was predicted that for a superconducting gap with BCS temperature-dependence Ao(r), the Josephson 
current through the TFT junction can reverse sign as the temperature T is raised. This unconventional temperature- 
dependence of I J is shown here in Fig. ([s]) . The reversal of Ij is a consequence of the splitting of the Andreev states 
by the transverse magnetic scattering terms in the barrier, which gives each state two distinct branches at ±|£^a| (ia) 
and ±|i?b| (±&). At T = 0, only the —a and ~b branches are occupied. Since \dEa/d(j)\ < \dEi,/d(f>\, the —6 branch 
makes the dominant contribution to the current and hence determines the direction of current flow. At any finite 
temperature the +a and -1-5 branches have a non-zero population, with the occupation of the -1-6 branch always larger 
than that of the +a branch since \Ea\ > \Ei,\. As the derivative of the +a and +b branches are equal and opposite to 
that of the —a and —b branches, respectively, this causes a reduction of the contributions to the current from both 
the a and b states. Because we have \Ea\ > \Ei,\, however, the change in occupation (and hence also the reduction 
in current) is greatest for the ±5 branches. If this reduction is sufficiently large the contribution of the a state may 
eventually dominate the current. If the sign of dEa/d(j> is opposite to dEi,/d(j), the direction of current flow then 
reverses. Since the a and b states are degenerate in the limits g —+ or a etc, it is thus not surprising that, for 
example, when = 7r/2 and 9 = g' = Z = the (i) sign change of Ij only occurs for sufficiently large g ^ I when 
a = [Fig. ([8])(a)], and (ii) is entirely absent for 7r/4 < a < 37r/4 [Fig. ^{h)]. 
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As the 0-dependence of the Andreev states is crucial for the temperature-dependent reversal of /j, it is of interest 
to consider the robustness of this effect in the presence of non-zero g' or Z. As shown in Fig. (|8|(c), both these 
terms suppress the reversal of /j, raising it to higher temperatures and reducing the maximum magnitude of the 
reversed current. The reason for this is due to the change in the 0-dependence of the Andreev states, as illustrated 
in Fig. (|8])(d). For both 5' ^ and Z 7^ the value of dEi,/d4> &t 4> — 7r/2 changes very little from the g' = Z — Q 
case; in contrast, the value of dEa/d<j> is considerably reduced. That is, the a state has a proportionately lower 
contribution to the current. To achieve a sign-reversal in Ij, therefore, a much higher temperature is required to 
sufficiently deplete the —b branch. More generally, the tendency of a finite Z or g' to enhance the T — current 
will suppress the sign-reversal effect: for example, when g = 2, a = O.IStt, and 9 ^ g' ^ 0, for Z > 1.5 we find 
that the sign of dEa/dcf) and dEb/dcj) are the same a,t cj) = tt/2, hence making a temperature-dependent reversal of Ij 
impossible (not shown). 

The two key requirements for a sign change in Ij with increasing temperature are that the Andreev states be 
non-degenerate and that the derivative with respect to (f> of the two distinct states have opposite sign. As such, when 
9 = this effect only occurs when the magnetic moment of the barrier has a finite transverse component. Because 
of the very general condition for its appearance, it is nevertheless not unique to the TFT junction considered here: 
a reversal of Ij with increasing T has also been predicted for a junction constructed by placing a magnetic barrier 
between two s-wave superconductors.'^ 



B. Non-aligned d- vectors 



We now consider the case when the d-vectors of the left and right superconductors are not aligned, i.e. 6* 7^ 0. This 



creates an interesting situation, as from Eq. (35) and Eq. ( |36| we see that the effective phase difference between the 
left- and right-superconductors is spin-dependent and given by — a9W^ We therefore first consider the case when 



g = 0, where the Andreev states are spin-polarized. The expression Eq. (60) for Ea thus becomes 



E, = kpAo^cos{[<l> - a9]/2), (62) 



with the same Dcr as in Eq. ( 60 1 . We sec that the effect of 6* 7^ is to shift the Andreev energies of each spin sector 
by 29 relative to one another [see black solid line in Fig. (|9|(a)]. This produces two sets of zero-crossings, at 

(j)n,± = {2n-l)TT±9, neZ, 

which are both evidenced in the current by jump discontinuities [see Fig. (|9|(d)]. Including now also a finite Z, the 
spin-t and spin-| Andreev states have different amplitudes [red broken line in Fig. (|9])(a) and Fig. ([9])(c)]. As in the 
case for 9 = 0, the Josephson current is spin-polarized. Here, however, it is possible that a spin current flows even 



when the charge current is vanishing. In such a case it is not possible to use the concept of the polarization [Eq. (61 )]; 
rather, we calculate the z-component of the spin current!^ 

Ijz = -^^iIjT^Ijl)- (63) 

For the states shown in Fig. ([9]) (a], we plot the corresponding Ijz in Fig. ([9|(b). Comparing the spin current with 
the charge current shown in Fig. ([9|(d), we see that for g' = 1, Z = g = 0, a, spin current flows a,t (f> = rnr even when 
Ij is vanishing. We find discontinuous jumps in the spin current occur at the zero-crossings of the Andreev states, as 
the sign of the Josephson current through the spin-cr states reverses as (f) is increased past (j)n.a — (2?! — 1)tt + (j9. In 
the case Z = 0, the sign of the spin current changes across the discontinuity; for Z = 1 the spin current is negative 
on both sides of the jump discontinuity at 4'n,+- The discontinuous sign changes in Ijz as a function of 4> may be 
thought of a "spin switch" effect, in analogy to the current switch effect discussed below. 

Because of the non-alignment of the d-vectors in the two superconductors, we see that if Z, 5' 7^ a current flows 



even when </> = 0, as clearly evidenced by Eqs. (57) and (58). Although it is still present when the barrier has a 
transverse magnetization, this result can be most easily understood when g = 0: in the spin-| sector there is an 
effective phase difference of —6* and an effective barrier oi Z — g'; in the spin-| sector there is an effective phase 
difference of 9 and a barrier Z + g' . Since the effective phases are of opposite sign, the Josephson current also flows 
in opposite directions in each spin sector. Assuming Z,g' > we see that as the effective barrier in the spin-| sector 
is larger than in the spin-f sector, the magnitude of the Josephson current through the former is smaller than that 
through the latter and hence there is a net current flow. This explanation is not readily applicable to the .9 7^ 
case as then the z-component of spin is no longer a good quantum number, but we nevertheless speculate that the 
a and h states retain some of the character of the g = states at finite g. As g is increased the mixing of the two 
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(e) '^l^ (f) 



FIG. 9: (color online) Comparison of the ^-dependence of (a) the Andreev states and (b) the corresponding Josephson spin 
current both with and without a finite Z for g' = 1, g = Q and 6 — n/3. In (c) we show the effect of a finite g on the g' = Z = 1 
Andreev states. We take a = it/4: and 9 — tt/S. (d) The Josephson current flowing through the Andreev states shown in (a) 
and (c). ^-dependence of (e) the Andreev states and (f) the corresponding Josephson current for g = 1.0, 9 = vr/S, a = 7r/4 
and at most one other barrier term non-zero. 



spin sectors also increases and the current at = is suppressed. The use of misaligned d-vectors to produce a 
supercurrent is similar in spirit to Asano's proposal to produce a spin current only at = OpS which is also present 
here [see Fig. (|9])(b)]. A charge current at 4> — was predicted for a junction between a p-wave and an s-wave 
superconductor with spin-orbit coupling in the tunneling barrier although our proposal involves a very different 
junction geometry, for the charge current to flow it similarly requires the breaking of time-reversal symmetry by the 
spin-dependent effective phases and barriers. 

Another remarkable aspect of the results for 9 ^ nn and Z, g' ^ is that the Andreev states are no longer symmetric 
about (j) = nn. As shown in Fig. (|9])(d), this has the striking result of removing the anti-symmetry of Ij about = nn, 
which is present in all other cases. This asymmetry is also found in the presence of a finite g (i.e. all barrier parameters 
non-zero). Despite significant changes in the 0-depcndcncc of the Andreev states in this case [blue dot-dashed line 
in Fig. ([9|(c)], the zero-crossings at 0ri,± remain: as is demonstrated by Eq. (50), the condition for zero energy states 
is independent of the barrier. 

We now turn to the case when g ^ 0, and at most one other barrier parameter is non-zero. Comparing Fig. ([9])(e) 
with Fig. ([2])(c), we see that in addition to producing zero crossings at 4>n,±, a finite 9 also removes the level crossings 
at = (j)LC when a ^ ac] for a — ac it also removes the level degeneracy (not shown). The effect of .Z 7^ or g' 7^ 
is qualitatively similar to that for 9 — 0: for both we observe a change in the sign of the curvature of Ea around 
(j) — (j)zci E^iid for Z ^ there is also a widening of the gap at (jj — (f>LCi while for the case 5' 7^ there is a substantial 
narrowing of the range of allowed energies. Because the zero crossings of the Andreev states are controlled only by 
the relative orientation of the d-vectors, jump discontinuities in Ij are always present when 9 ^ nn [see Fig. (|9])(d) 
and Fig. ([9])(f)]. As was discussed in Sec. IIIAl for 9 — 0, when the a and b states are non-degenerate, the b state 



in general gives a greater contribution to Ij than the a state; numerical investigations show that this is also the case 
for 6^0. Because the contribution to the current from the b states reverses sign for 0„^_ < < 4'n,+ , the jump 
discontinuities in Fig. ([9])(f) are therefore accompanied by a reversal of the sign of Ij. Furthermore, we note that the 
dependence of Ij on Z and g' for 4>n,- < 4> < 4'n,+ is very different to that for < (j) < </>„._ : for the latter 

we see an enhancement of the current, whereas |//| is reduced in the former case. This is due to the changes in the 
0-dependence of the Andreev states, which tend to enhance the contribution to Ij from the a state in the region 
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FIG. 10: (color online) Comparison of (a) 7j and (b) the Andreev states for Z = and Z — Zci = 1.732051 . . ., with fixed 
9 — tt/3, g — 2.0, g' = and a — ■k/2. The dotted line in (a) is a guide to the eye. The current at fixed is plotted as a function 
of Z in (c) for the same parameter values as before. Note the sign-change at the two critical values Zd and Zc2 = 2.88675 . . .. 



(j^n,- < </> < (/'n.+ i while leaving the contribution from the b state mostly unaffected. As the contribution from the a 
and b states is of opposite sign for these values of 4>, this leads to an overall decrease in \Ij\. The current increase for 
. is due to the subtle effects discussed in Sec. 



Ill A 



(Pn-1,+ <</)<< 

A remarkable example of the competing contributions from the a and b states is found in Fig. (10): for given 9, 
it is sometimes possible to select barrier parameters such that the current from the a and b states completely cancel 
each other for a finite range of ( 
condition 



Analysis of Eq. (58 1 reveals that this only occurs for g,Z ^ and g' ~ 0, when the 



D [g^ cos{9 - 2a) - (1 + Z^) cos 9] - sgnS = 



(64) 

is satisfied. As presented in Fig. ([10|(a) for 9 — 7r/3, g — 2.0 and a ~ 7r/2, a barrier potential of approximately 
Z = Zci = 1.732051 . . . leads to a cancellation of the current for (\)n- < <j) < <j>n,+- The current cancellation refiects a 
special relationship between the a and b states at Zd'- for this choice of barrier potentials the moduli of the energies 
of the two states are equal up to a ^-independent constant C, i.e. 



< (p < 

< (I) < (Pn,+ 



(65) 



Note from the definition Eq. (55) that Ea{(t>) is always positive. When the contribution to the current from the b 
state reverses sign for 0„__ < <j) < 4'ri,+ : it exactly cancels the contribution to the current due to the a state. Further 
increasing Z above Zd , the contribution to the current from the a state dominates that from the b state and hence the 
sign of the current for (pn,- < 4> < 4'n.+ is reversed compared to the Z = case, as shown in Fig. (lOl(c). Although 
the current still displays discontinuities at (f>n.±, it does not change sign across the intervals ((/>n,± — 0~^-,4'n,± + 0+). 
As Z is increased a second zero current state is found at Z — Z^i = 2.88675 . . ., above which the b state again has 
the dominant contribution to Jj and the sign of the current for (\)n.- < (t> < 4>n,+ is the same as at Z = 0. We note 
that since lim^^oo Ij = 0, the magnitude of the current for (jin,- < (j) < (t>n,+ will take a maximum at some Z > Zc2, 
beyond which it asymptotically decreases (not shown). 



1. Rotating the d-vectors: Equilibrium limit 

The jump discontinuity in Ij observed in Fig. (|9]) can be used to construct another Josephson current switch. For 
fixed phase difference between the left and right superconductors, rotating d/j will cause a zero-crossing of the a states 
at 

9n,± = {2n+l)7r±<f>, n e Z, 



and hence a jump discontinuity in the current as presented in Fig. ( 11 ). This clearly forms the basis of a switch effect, 
as for 9 close to ^„.±, small changes between the relative alignment of the d-vectors can switch the junction between 
different current states with opposite direction of Ij. The robustness of this switch effect, however, depends upon 
the ratio of the period of the rotation Tg to the relaxation time r. In this section we consider the equilibrium limit 
where Tq ^ t and so the system is always in thermal equilibrium; the opposite limit where Tg <^ t is considered in 
the following section. In the equilibrium limit, the occupation of the Andreev levels is given by the Fermi distribution 
throughout the period of the rotation. 
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FIG. 11: (color online) Current as a function of 9 for various barrier parameter combinations with g = 1.0 at </!> = 7r/5 and (a) 
a = 0; (b) a = 7r/3; (c) a = 7r/2. (d) Current as a function of for various barrier parameter combinations with g' — 1.0 at 
(j) — 7r/5 and a — 0. 



In Ref. 122] it was noted that the 



dependence of the two current states can be significantly altered by changing the 
value of a [presented here in Fig. ( 11 1 (a) through (c)]. In particular, note that only at a = and a = Ucis the current 
symmetric with respect to about 9 — nn; for a intermediate between these values, the current is strongly skewed, 
varying almost linearly with 9 in Fig. pT])(b). The ^-dependence of the current can also be modified by a finite Z or 



g', as seen for example in Fig. (11 1 (a) where we observe that the curvature with respect to 9 of the "positive" and 



"negative" current states can reverse sign as Z or is increased. As a is increased towards a^, however, the changes 
produced by non-zero Z or g' become much less pronounced [e.g. Fig. ( 11 l(c)]. As shown in Fig. ( 11 )(d), more radical 
alterations in the 0-dependence of Ij can be obtained when g' and Z are both finite. 

Close inspection of Fig. ([TT])(b) for Z — 1 shows that the current changes sign at 9 ^ O.ln, slightly before the jump 
discontinuity aX 9 — O.Stt; also when Z — \, a sign change is observed in Fig. (ll)(d) at 9 m O.Stt both with and 
without a finite g. These results indicate that the switch effect can be spoiled by sufficiently large Z: although a 
jump discontinuity between the two different current states still o ccur s at 0n,_, the current does not reverse sign as we 
move from 9n- — 0^ to 6'„__ + 0+. The similar behaviour in Fig. ( 11 l(b) and (d) nevertheless arise from qualitatively 
different ^-dependencies of //, which can be seen in Fig. (12 1. We first consider the current corresponding to the 
Z = 1 case in Fig. (ll)(b). Since g' = 0, Ij is antisymmetric about = mr: as such, the continuous reversal of Ij 
at (f> = 7r/5 with increasing 9 implies that there must be an intermediate state where Ij = for < (j) < (pn,- 

[the reverse of the situation shown in Fig. (10 1 (a), where /,/ = for (j)n- < (j) < This is indeed observed 

in Fig. (12 1 (a), where the current for (pn~i.+ < (f> < (f>n,- vanishes at ~ 0.7157r. When g' and Z are both finite, 
however, Ij is not antisymmetric about (p — nir as shown in Fig. (12)(b). The sign change in Fig. (lll(d) with 
increasing 9 is therefore accomplished by a downward shift in the value of the current for <j)n~i,+ < (j) < 4'n,-'- for 
0.6487r < 9 < TT, the current is only negative for this range of cj) values. Despite these differences, the presence of a 
finite Z is apparently a prerequisite: exhaustive numerical investigations for Z ~ have failed to find an example 
where the sign of Ij does not reverse across the discontinuity. 



2. Rotating the d-vectors: Adiabattc linut 



The adiabatic limit is reached when the rotation period of dn is much faster than the relaxation time of the system, 
i.e. Tg <^ T. As was argued in Ref. |22l the occupation of the Andreev states do not in this case assume their 
equilibrium values throughout the entire rotation, but rather remain fixed at their initial levels. This implies that 
the discontinuous jumps in current found in the equilibrium limit due to the zero-crossings are absent: as displayed 
in Fig. (13 1 (a), starting from an initial state at 9 = 0, the —a level remains fully occupied for all 9, even when it 
is higher in energy than the +a level for 6'„_ < 9 < 9n+- Consequently, the current is a continuous function of 9 
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FIG. 12: (color online) Current as a function of <j) for Z = 1, a — n/3 for various values of 6, and for (a) g — 1.0, g' = 0, 
a = ^/3; (b) g = 0.0, g' = 1.0. 




FIG. 13: (color online) (a) Occupied and unoccupied Andreev states during the adiabatic rotation of d_R. We g — 1.0, 
g' = Z = 0, (j) = n/5 and a = 0. (b) Josephson current as a function of 6 for rotations in the equilibrium and adiabatic limits. 
Parameters are the same as in (a). 



and there is no change in the sign of // during the rotation. This is compared to the result in the equilibrium limit 
in Fig. ([13]) (b). 

A reversal of the Josephson current as dj^ is rotated is still possible in the adiabatic limit, although this requires a 
finite potential barrier Z and is not related to the zero energy crossings of the Andreev states. As seen in Fig. ( 11 l(b) 
for Z = l.O and a — 7r/4, the Josephson current changes sign before the discontinuous jump is reached. As such, 
starting at 6* = and rotating adiabatically we therefore encounter a reversal of Ij sX 9 k, O.Gtt. The current 
changes back to its initial direction at 9 k, I.Itt (not shown). 



C. Critical Current 



The critical current Ijc is defined as the maximum magnitude of the current that can be carried by the Josephson 
junction. The maximum with respect to 4> can occur at three places: a stationary point of /j, or on either side of the 
jump discontinuity at 4>n,± (i-e. at (j)n,+ i 0"*" or equivalently ^pO"*"). These three critical current locations may be 
thought of as defining "phases" of the junction. As the barrier parameters or the relative orientation between the two 
d- vectors are varied, we find "transitions" between the different phases. An example of this is shown in Fig. ([9])(d): 
as Z or g' is increased from zero, the position of Ijc changes from + 0+ to the stationary point. The phase 
boundaries are marked by lines of non-analyticity in Ijc- For simplicity, below we specialize to the case g' — Q. 

In Fig. (14 1 (a) we plot Ijc in the Z-g plane for fixed 9 — O.SStt and a — Q.2tt. A line of first-order non-analyticity, 
gci{Z), is shown as the solid white line; the broken white line shows a line of third-order non-analyticity at gc2{Z) < 
gci{Z). The line gc2 defines the boundary between the phase where Ijc is located &t cf) — — 0^ (ff < 9c2), and 
the phase where the maximum is located at the stationary poin t {gc2 < 9 < ffci)- For g > gd, the current maximum 
is always located a,t cj) = (pn- +0^. As can be seen in Fig. (14l(b), the line of third-order non-analyticity gc2 is not 
easily detectable in a plot of Ijc as a function of g; as such, below we only consider lines of first-order non-analyticity. 

We plot the line of first-order non-analyticity in the Z-g plane for different values of 9 in Fig. ( 14 1 (c) . For orthogonal 
d-vectors (i.e. 9 — tt/2) and a — tt/5, we do not find any non-analyticity in the current (since only the squares of Z 
and g enter into Eq. ( 56 ) , we need only consider the case Z, g > 0) . Upon an infinitesimal decrease in the value of 
9, a line of first-order non-analyticity appears at 5 = 0+; further decreasing 9, we find that this line moves to higher 
values of g. As 9 0, the line moves to infinity. In contrast, for Z,g < 2 the Z-g plane remains free of first-order 
non-analyticities as 9 is increased from tt/2 to tt: we have not found any first-order non-analyticities up to Z,g = 10, 
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FIG. 14: (color online) (a) Critical current as a function of Z and g for fixed 9 = O.SStt and a = 0.2tt. The value of Ijc/ekpAo 
is given by the scale on the right. The white lines denote boundaries between the different phases of the junction: Ijc occurs 
at (j>n,± ± 0^ in phase A, at the stationary point in phase B, and in phase C it occurs at (j)n,± =F 0^. The solid line indicates 
a first-order non-analyticity in Ijc; the broken line indicates a third-order non-analyticity. (b) Critical current as function of g 
for fixed Z. The point of first-order non-analyticity is clear as the sharp kink; the position of the third-order non-analyticity is 
indicated by the arrows. 6 and a are as in panel (a), (c) Lines of first-order non-analyticity of Ijc in the Z-g plane for a — n/5 
and various values of 9. (d) Lines of first-order non-analyticity of Ijc in the a-g plane for Z = 1.0 and various values of 9. 



beyond which it becomes difficult to numerically check for non-analytic behaviour. This curious dependence upon 9 
can be better appreciated by examining the line of first-order non-analyticity in the a-g plane [Fig. ( |14[ )(d)]. We find 
here an interesting symmetry about a = 0: our numerical investigations indicate that the line of non-analyticity for 
a di^-vector orientation tt/2 < < tt can be found by reflecting the line of non-analyticity for d/j- vector orientation 
TT — 6 about a — 0. In the limit 9 — {2n + 1)tt/2, the lines of non-analyticity are found at constant a — n7r/2. We 
have also numerically checked up to high accuracy that the curves in the a-g plane are symmetric about the line 
a = 6/2 {a = tt/2 — 6/2) for < 9 < tt/2 {tt/2 < 6 < tt). Note that the plot in the a-g plane is periodic in a with 
period tt. The constant line a — tt /h is drawn in Fig. ([T4|(d) as the thin dotted line; this corresponds to the situation 
in Fig. (14l(a-c). We see that as 9 is increased from zero to 7r/2, the intersection of the lines of non-analyticity with 
a = tt/5 occur at progressively smaller values of g, consistent with Fig. (14)(c). As we cross to the tt/2 < 9 < tt case, 
we see that in the range < g < 5 there is no intersection of the lines of non-analyticity with a — tt/5. Thus, 



the 

absence of non-analyticity in Fig. p4|(c) for tt/2 < 9 < tt can be explained by the a-dependence shown in Fig. ( 14)(d). 

In Fig. (15) (a) we plot Ijc as a function of 6 and g, with fixed Z = 0.5 and a — O.Stt. The line of first-order non- 
analyticity is again shown by the solid white line, while the third-order non-analyticities are indicated by the broken 
white lines. Note the different phases on either side of the first-order non-analyticity as compared to Fig. (14 1 (a). In 
general, a line of first-order non-analyticity can separate any combination of different phases. The critical current as 
a function of 9 is shown at several values of fixed g in Fig. (15)(b). Although the first-order non-analyticity is clearly 
evident as the abrupt change in the slope of Ijc, the third-order non-analyticities do not clearly coincide with any 
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FIG. 15: (color online) (a) Critical current as a function of Q and g for Z = 0.5 and q = O.Stt. The value of Ijc/ekp^o is 
given by the scale on the right. The white hues denote boundaries between the different phases of the junction: Ij^ occurs at 
0n,± ± 0^ in phase A, at the stationary point in phase B, and in phase C it occurs at 4'n,± =F 0^. The solid line indicates a 
first-order non-analyticity in the broken lines indicate third-order non-analyticities. (b) Critical current as function of 6 
for fixed g. The point of first-order non-analyticity is clear as the sharp kink; the third-order non-analyticities are indicated by 
the arrows. Z and a are as in panel (a), (c) Lines of first-order non-analyticity of Ijc in the 6-g plane for Z — 0.5 and various 
values of a. (d) Lines of first-order non-analyticity of Ijc in the 9-Z plane for g = 2.0 and various values of a. 



noticeable feature. 

The line of first-order non-analyticity in the 9-g plane displays a complicated variation with a as shown 
in Fig. (15l(c). We see that the limit a = mr is somewhat pathological, as here we have a line of non-analyticity at 
6 = 7r/2 which branches into three lines at 5 w 1.2, with the top and bottom branches respectively converging to 
9 — TT and 9 — as g ^ 00. For finite a ^ a^, the intersection of the branches disappears and the lines of first-order 
non-analyticity for tt/2 < 9 < tt and < 9 < tt/2 are unconnected. The it /2 < 9 < n branch disappears as a 
is increased, while the g 00 limit of the bottom branch moves up to = 7r/2. At a — ac, a line of first-order 
non-analyticity is found only at constant 9 — Tr/2. Interestingly, this is also the line of non-analyticity when transverse 
magnetic terms are absent in the barrier. This can be seen by examining the lines of non-analyticity in the 9-Z plane 
for fixed g = 2.0, as shown in Fig. ([T5])(d). In the limit Z/g — > 00, the potential terms in the barrier Hamiltonian 
dominate the transverse magnetic terms: we accordingly find that for any value of a, the < 9 < 'k/2 branch of 
non-analyticity asymptotes to 6* = tt/2 as Z ^ 00. The 7r/2 < 9 < tt branch, in contrast, is always suppressed beyond 
some finite Z. 
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IV. CONCLUSIONS 

In this paper we have have presented an extended analysis of the novel TFT junction first considered in Ref. [22] 
Our aim has been two-fold: by including a potential scattering term or a component of the magnetization normal to 
the barrier interface, we have attempted to assess the persistence of the many new effects identified in Ref. under 
a more general description of the system; we have also investigated the possibility that the presence of these extra 
terms in the barrier Hamiltonian gives rise to the emergence of other unconventional Josephson behaviour, absent in 
the previously-studied case. 

We predict that the additional degrees of freedom in the description of the barrier cause three new effects at T = 0: 
the spontaneous generation of a Josephson current by misalignment of the d-vectors, even when there is no phase 
difference between the condensates; the existence of a special line in the parameter space such that Ij vanishes over 
a finite range of 0; and that the current through a magnetic barrier with potential scattering can show a strong 
enhancement over its value in the absence of a potential term, no matter what the orientation of the ferromagnetic 
moment is. The z-component of the spin current across the junction was calculated for the case g = Q, and it was 
found that it could be finite even when the charge current is vanishing. Furthermore, the critical current through 
the TFT junction was studied as a function of the barrier parameters and d- vector alignment. The location of the 
critical current along the Ij vs. (/) curves was classified into three qualitatively distinct categories. This enabled the 
construction of "phase diagrams" for the TFT junction, where the different "phases" correspond to different locations 
of the current maximum. In certain cases a first-order non-analyticity in the critical current is found at the phase 
boundaries; we propose using these non-analyticities as a test of our theoretical knowledge of the TFT junction. 
Our results also indicate that the current switches identified in Ref. 22 are largely robust to these alterations in the 
properties of the barrier. Even when g' and Z are comparable in magnitude to g, we find that it is still possible to use 
the transverse component of the barrier magnetization to tune the system between "off" (/j ~ 0) and "on" (/j ^ 0) 
current states, or to create abrupt reversals in the current direction by small changes in the relative orientation of 
the d-vectors. On the other hand, because of the sensitive dependence of Ij upon the ^-dependence of the Andreev 
states, the reversal of the current with increasing temperature^ is strongly suppressed by a finite Z or g' . 

A natural question concerns the observability of the effects predicted above. We first note that our analysis is 
most applicable to systems where the BdG description of the quasiparticle excitations in the superconducting state is 
expected to be good. This should be the case for Sr2Ru04, where the superconducting state is well-described by the 
Balian and Werthamer generalization of the BCS theory to triplet pairing.f^lSlfil The application to the p-wave state in 
heavy-fermion or organic superconductors is less certain, as it is not c lear t o what extent the standard weak-coupling 
description of the quasiparticle states is reasonable in these materials.l^SlIl] 

More difficult is the question of experimental control over the different junction parameters. Although it is impossible 
to forecast future developments in the fabrication of quantum devices, we expect that the orientation of the barrier 
magnetic moment and the phase difference (j) should be the easiest parameters to manipulate. Much more challenging 
would be the orientation of the d-vectors: as these point in a fixed direction in the crystal, the experimentally- 
accessible values of 9 would be limited by the different orientations that one could grow a Sr2Ru04 crystal upon the 
substrate provided by the magnetic barrier. The increasing degree of control in growing oriented crystal interfaces 
for tunneling experiments in cuprate superconductors makes us hopeful that this is not an insurmountable obstacle 
to the experimental study of the ^-dependence of the current .1^^ Alternatively, one could imagine a purely mechanical 
control over one of the superconducting slabs: indeed, there has recently been much interest in integrating mechanical 
degrees of freedom into superconducting tunneling experiments, albeit thus far in terms of "charge shuttle" effects.!^ 
In any case, we emphasize that the experimental realization of the device proposed here would provide important 
insights into the interdependence of ferromagnetism and triplet superconductivity. 
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APPENDIX A: BOUNDARY CONDITIONS 



Here we provide explicit expressions for the component matrices of A/tft in Eq. (48 1. Assuming that Pi^ ~ k,y ^ Ki,, 
we obtain 
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In Eqs. (A3l and (A4|, r = k^/kj^ is the ratio of the Fermi vectors in the left- and right-superconductors. 



APPENDIX B: TFT JOSEPHSON JUNCTION FOR REDUCED SCATTERING TERMS 



Here we specialize the TFT junction solution to Z = g' = 0, i.e. the barrier is composed only of a transverse 
magnetic scattering term. We also assume that r — 1, Al — Afi = Aq and niL = mn = m. This is the case treated 
in Ref. [22] and here we intend to make contact to the slightly different form of the resulting Andreev energies given 
there. We show below analytically that in this case always both solutions with Andreev state energies Ea^b exist (we 
note in passing that a numerical investigation into the case where also g' and Z are allowed to be non-zero shows that 
the same is likely true although we have not succeeded in proving this analytically). 

Eq. ( 49 1 now becomes 



(e2*'^/64)detMTFT = cos^ 27 - sin^ - 2cos 27COS 0cos(6'l - 9r) + cos'^{0l - Or) 
—Ag^ cos^ 7 [cos^ 9l+ cos^ Or — cos 27 — cos cj) cos{Ol + Or)] 
+4:g* cos'* 7 = 0, 

with Ol and Or as defined after Eq. (49). As before, this may be recast as a quadratic equation for E^: 
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The non- negative solutions of Eq. ( B2 1 are 



Ea{b) = kpAoVD 



(B7) 



which is easily seen to coincide with Eq. (4) in Ref. [22] , 

In this case, we may prove that the solutions exist for all parameter values. For this to be true, we require that 



Since 



and 



(i) : A-C>0, 

(ii) : B-C>0. 
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(A - C)\^^^ = cos'{0L - Or) cos' \^ > 0, 



(Bll) 



and since < £> < 1 and since A — C is linear in D, we have A — C > for all admissible D, i.e., for general g and 
thus (i) is always true. Since 



(B - = sm'{0L - On) sm\dL + Or) + sin^ \(j)Cos\6L + 0/?) > 



and 



{B - = sin^OL - Or) sin^ i0 > 



(B12) 



(B13) 



and since < I? < 1 and since S — C is linear in D, we have B — C >Q iox all admissible D, i.e., for general g and 
thus (ii) is always true. Since (i) and (ii) always hold, Andreev states with energies Ea^b always form in this case. 
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